Determining the current polarization in Al/Co nanostructured point contacts 
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We present a study of the Andreev reflections in superconductor/ferromagnet nanostructured 
point contacts. The experimental data are analyzed in the frame of a model with two spin-dependent 
transmission coefficients for the majority and minority charge carriers in the ferromagnet. This 
model consistently describes the whole set of conductance measurements as a function of voltage, 
temperature, and magnetic field. The ensemble of our results shows that the degree of spin polar- 
ization of the current can be unambiguously determined using Andreev physics. 

PACS numbers: 74.45.+C, 72.25.-b, 74.78.Na 



The field of spintronics is largely based on the abil- 
ity of ferromagnetic materials to conduct spin-polarized 
currents 0. Thus, the experimental determination of 
the degree of current polarization has become a key is- 
sue. Recently the analysis of Andreev reflections in su- 
perconductor/ferromagnet (S/F) point contacts has been 
used to extract this spin polarization in a great vari- 
ety of materials 0, 0, 0> IE IE 01 • The underlying idea 
is the sensitivity of the Andreev process to the spin of 
the carriers, which in a spin-polarized situation is man- 
ifested in a reduction of its probability [8j. The theo- 
retical analysis of these S/F point-contact experiments 
has been mainly carried out following the ideas of the 
Blonder-Tinkham-Klapwijk (BTK) theory 0. Different 
generalizations of this model to spin-polarized systems 
have been proposed, in which with an additional phc- 
nomenological parameter P, the spin polarization of the 
ferromagnet, excellent fits to the experimental data have 
been obtained 0, & IE IE • However, a microscopic 
justification of these models is lacking [1(1 lljl U2 • ■^• e " 
cently, Xia el al. 0] have combined ab initio methods 
with the scattering formalism to analyze the Andreev 
reflection in spin-polarized systems. Their main conclu- 
sion is that, in spite of the success in fitting the experi- 
ments, these modified BTK models do not correctly de- 
scribe the transport through S/F interfaces. Therefore, 
at this stage several basic questions arise: what is the 
minimal model that describes on a microscopic footing 
the Andreev reflection in spin-polarized systems? And, 
more importantly, can the current polarization be exper- 
imentally determined using Andreev physics? 

In this paper we address these questions both experi- 
mentally and theoretically. We present measurements of 
the differential resistance of nanostructured Al/Co point 
contacts as a function of voltage, temperature, and mag- 
netic field. To analyze the experimental data we have 
developed a model based on quasiclassical Green func- 
tions, the main ingredients of which are two transmis- 
sion coefficients accounting for the majority and minor- 



ity spin bands in the ferromagnet. We show that this 
model consistently describes the whole set of data, which 
unambiguously demonstrates that the spin polarization 
of current in a ferromagnet can indeed be determined 
employing Andreev reflection. 

We have fabricated Al/Co point contacts following 
the process described in Ref. 14]. Briefly, a bowl- 
shaped hole is drilled through a 50-nm thick silicon ni- 
tride (Si3+ a; N4_ 2 ;) membrane by means of electron-beam 
lithography and reactive ion etching. The smallest open- 
ing in the insulating membrane has typically a diameter 
of 5nm. Finally, 200 nm of Al and d Co = 6, 12, 24 or 
50 nm of Co plus (200 nm - dc ) of Cu are deposited 
by electron-beam evaporation under ultra-high vacuum 
conditions (~ 10~ 9 mbar) on each side of the membrane. 
A schematic of the samples is shown in Fig. ^a). The 
differential resistance R was measured with lock-in tech- 
nique in a dilution refrigerator. A dc current was super- 
imposed to the small measuring ac component and both 
R and the voltage drop V were recorded simultaneously. 

As a reference we show in Fig. ^b) the Andreev spec- 
trum, i.e. the differential conductance G as a function of 
the voltage V, of an Al/Cu sample. In all the spectra in 
this paper, G and V have been normalized by the nor- 




Co Cu eV/A 

FIG. 1: (a) Schematic of an Al/Co nanocontact. (b) Andreev 
spectrum of an Al/Cu contact at 95 mK (black circles). The 
dashed line is the fit obtained with the BTK theory [jj yield- 
ing the transmission r = 0.781 and the gap A = 206 fieV. 
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12 
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199 
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12 


13.3 


100 


188 
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0.970 


0.40 


#5 


24 


6.00 


98 


180 


0.382 


0.989 


0.44 


#6 


24 


3.58 


97 


193 


0.399 


0.983 


0.42 


#7 


50 


15.7 


99 


172 


0.370 


0.994 


0.46 


#8 


50 


3.59 


97 


198 


0.392 


0.986 


0.43 



TABLE I: Transmissions, Tx,i, polarization, P, and gap, A, 
for the Al/Co samples as determined by a fit of the Andreev 
spectra for T ~ 100 mK with our model. 



FIG. 2: Andreev spectra of four Al/Co point contacts with 
different Co film thickness dco- The solid line is a fit to the 
data with our model (see Table [IJ. 



mal state conductance Gn and by the zero-temperature 
superconducting gap A of the Al electrode, respectively. 
Gjv showed to be completely independent of V in the 
range eV < 5—10 A. Since the estimated mean free paths 
of the Cu and Al electrodes are ~ 60 nm or longer at low 
temperatures, all the contacts studied are in the ballistic 
regime. In the Al/Cu case (Fig. H^b)) the BTK theory 
fits the experimental data very well (see figure caption for 
details). In the case of Al/Co, the ferromagnetic layer 
causes a reduction of the Andreev spectrum amplitude 
as compared to the Al/Cu contacts (see Fig. 01. Notice 
that, although both the normal state resistances and the 
Co layer thicknesses of the samples differ strongly (see 
Table [IJ, the Andreev spectra are all quite similar. This 
indicates that we are observing an intrinsic property of 
Al/Co point contacts. 

The minimal model necessary to describe transport 
in S/F contacts should account for the spin-dependent 
transmission, which is inherent to any junction where fer- 
romagnets are involved. We have developed a model that 
fulfills this requisite in the framework of the quasiclassi- 
cal Usadel theory [lflL flfil ] , describing a system in terms of 
two retarded Green functions, g(r,e) and f(f, e), which 
depend on both space and energy and satisfy g 2 + f 2 = 1 . 
For transport through interfaces this theory must be sup- 
plemented with boundary conditions, which can be for- 
mulated in terms of a normal-state scattering matrix S. 
Our choice to model an S/F interface is given by (we 
restrict ourselves to a single conduction channel) 



S = 



t = 











o 



(i) 



where t^i and r^j, are the spin-dependent transmission 
and reflection amplitudes, respectively. The transmission 



coefficients r-f^ = \t-\.i\ 2 are the central quantities of our 
model. They contain the microscopic properties relevant 
for transport, i.e. the spin-split band structure of the fer- 
romagnet, the electronic structure of the superconductor 
and the interface properties. 

The current Isf through the S/F point contact is com- 
puted following standard procedures 0|- ^ can be sep- 
arated in two spin contributions, Isf = Ij + /j., where 
each can be written in the BTK form 9) 



la = 



de [n F (e - eV) - n F {e)} [1 + A a {e) - B a {e)} . 

(2) 

where np is the Fermi function, and A a (e) and B a (e) are 
the spin-dependent Andreev reflection and normal reflec- 
tion probabilities, respectively. These are given by A„ = 
r a T^\f/V\ 2 and B a = |(r CT + r_ ff ) + (r a - r_ a )g\ 2 /\V\ 2 , 
where r CT = ^Jl ~ T a an d V = (1 + r a r- a ) + (1 - r a r- a )g. 
The Green functions are evaluated right at the interface 
at the superconducting side. In the point-contact geome- 
try we can ignore the proximity effect, which means that 
g and / only contain properties of the superconducting 
electrode. In the case of a BCS superconductor in zero 
magnetic field g = ie/yj A 2 — e 2 and / = —i(A/e)g, and 
the zero-temperature conductance adopts the form |l8j 



^SF — —r- 



(i+nn) 

r T r i + ( r T+ r l~ 



-4r T r i ( eV/A) 2 

nn)Vi-(A/eV) 



eV < A 
eV > A 



(3) 

In the absence of spin polarization (tj = r| ) this formula 
reduces to the BTK result |9j . The normal state conduc- 
tance is given by Gn — (e 2 /h){T\ + rj.), and the current 
polarization is defined by P = |tj — tJ/(t| + t\). The 
main approximation of this model is the assumption that 
we can describe the point contact with a single pair of 
transmission coefficients, r^j, which will be finally justi- 
fied by the agreement with the experiment. 
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FIG. 3: (a) Andreev spectrum for sample #2 for different temperatures. For clarity, the curves are shifted downwards 
successively by 0.05 units with increasing temperature. The solid lines are the calculated spectra with our model, (b) 
Normalized resistance R/Rn as a function of temperature for three Al/Co samples. T is normalized to the gap A as obtained 
from the Andreev spectra (see Table 0. The curves are shifted downwards successively by 0.12 units. The red lines are the 
calculated R{T). For sample #2, the dashed line has been calculated including the effect of a residual magnetic field of 5 mT. 
As a reference, we also show this data for the Al/Cu contact of Fig.Q^b) (the theoretical result corresponds to the nonmagnetic 
BTK theory). The shaded region is covered by a set of curves given by {tj ±0.01,t; ±0.01} for sample #3. (c) Andreev 
spectrum for sample #4 measured at 100 mK for different magnetic fields. The curves are shifted upwards successively by 0.2 
units. The inset shows the zero-bias conductance as a function of the field. The critical field of the sample is (j,qH c = 15.0 mT. 
The red lines are the calculations using d/Ao = 3.8. 



As we show in Fig. El using tj^ and A as free parame- 
ters our model yields an excellent fit to the Andreev spec- 
tra of the Al/Co contacts for temperatures T ss 100 mK. 
These parameters for a total of eight contacts are listed 
in Table Their deviations from sample to sample are 
remarkably small, leading to small uncertainties in the 
mean values given by f T = 0.40 ± 0.02, rj. = 0.98 ± 0.01 
and A = (190 ± 10) /ieV. The total current is of course 
symmetric with respect to the exchange of r-j and tj,, 
which implies that we cannot assign a transmission co- 
efficient to the majority or minority charge carriers in 
Co. Nevertheless, we expect the high transmissive coeffi- 
cient rj, to correspond to the minority electrons, because 
of their higher density of states at the Fermi level corre- 
sponding to the Co 3d band. In our contacts the mean 
value of the current polarization is P = 0.42 ± 0.02. An 
analysis of our experimental data for T w 100 mK with 
the widely used model of Ref. gives fits of similar 
quality, but yields ~ 15% smaller values for P. It is im- 
portant to stress that this model cannot be mapped onto 
ours, it is not rigorously founded, and misses the funda- 
mental ingredient of a spin-dependent transmission. 

The rest of the paper is devoted to illustrate the con- 
sistency of the model, and in turn of the determination of 
the polarization P. We show that fixing the set {r|, rjj 
and A, as obtained from the spectra at T « 100 mK, 



the model describes without any additional fit parame- 
ter the temperature and magnetic-field dependence of the 
conductance. For instance, in Fig. |3fa) the temperature 
dependence of the Andreev spectrum of the sample #2 is 
depicted. As can be seen, the model describes the whole 
temperature range by simply using the BCS tempera- 
ture dependence of the gap. A more stringent test of our 
model is shown in Fig. |3Jb). Here, we compare the tem- 
perature dependence of the zero-bias resistance with the 
theoretical prediction. The agreement is excellent, apart 
from the deviations close to the critical temperature. We 
attribute them to the existence of a stray field (~ 5 mT) 
created by the Co film. This idea is supported by a calcu- 
lation (see below for details) of R(T) in the presence of an 
external field (Fig.|jjb)). It is worth stressing that R(T) 
is extremely sensitive to the transmission (see curve for 
sample #3 in Fig. EJb)), which illustrates the accuracy 
in the determination of {r-j-, t^}. 

We have also measured how a magnetic field, H, par- 
allel to the insulating layer modifies the Andreev spec- 
tra (see Fig.[2{c)). There are three main effects: (i) the 
height of the two maxima diminishes with increasing field 
and their positions are shifted to lower voltages, (ii) As 
can be seen in the inset of Fig. Etc), the zero-bias con- 
ductance is constant for fields below the critical field, 
(iii) The transition to the normal state is abrupt. To 
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understand these features we now study how the order 
parameter, A, is modified by the field. We use two ap- 
proximations: (a) in the Al electrode the mean free path 
(/ ~ 60 nm) is much smaller that the superconducting 
coherence length (£o ~ 300 nm), which justifies the use 
of the diffusive approximation (I -C £o) and the Usadel 
theory, (b) For our Al films £o is greater than the elec- 
trode thickness, d, which means that we can assume that 
A and the Green functions are constant throughout the 
sample. With these approximations the Usadel equation 
reduces to the generic equation that describes the effect 
of different pair-breaking mechanisms such as magnetic 
impurities, supercurrents or magnetic fields 



e+iTg{e,H) = iA 



where T 



2De 2 
he 2 



(A 2 ), (4) 



where D is the diffusion constant, T is a depairing en- 
ergy, which contains the effect of the magnetic field, 
and (A 2 ) is the average value of the square of the vec- 
tor potential along the thickness of the Al film. Ad- 
ditionally, the order parameter A must be determined 
self-consistently jj^- in Al the London penetration 
depth is typically Ao ~ 50 nm, which in our case is 
smaller than the thickness, d. This implies that the 
external field is partially screened inside the sample. 
Thus, the vector potential appearing in Eq. (4) must 
be determined solving the Maxwell equation: V 2 A — 
— (Ait/c)j, where j is the supercurrent density given by 
J(f) = -(2a N /Hc)A(r)J™ de tanh (/3e/2) lm(/ 2 ), where 
ctjv is the normal conductivity of the Al sample and 
/3 = (ksT)" 1 . The solution of the Maxwell equation 
yields the following expression for the depairing energy 



T(H) 



Go 



r 2 cosh (r/2) \ r 



sinh(r) 



De 2 d 2 H 
6hc 2 



(5) 

1/2 



where r = (d/X ) [(2/tt) /" de' tanh (&f\ Im(/ 2 ) 
Here, the prime indicates that the energy variables 
are measured in units of the zero-temperature gap in 
the absence of field, A , and A = \/hc 2 / (An 2 a n A ) . 
In Eq. (5) a is the pair-breaking parameter for a 
thin film [19j, which can also be written as a/Ao = 
(l/127r) [Hd/ HcbXo] , where H c b is the bulk critical field. 
For Al noHcb = 9.9 mT. Notice that the ratio d/Xo is 
the only parameter that enters our analysis. Since d/Xo 
determines the critical field of the Al films, H c , we fix its 
value by means of an independent measurement of R(B) 
at T « 100 mK. For our samples, we find H c w 1.5i? c h, 
which in our theory corresponds to d rs 4Ao- Thus, us- 
ing Eq. (2) with the self-consistent solution of Eq. (4) 
for the Green functions, we calcute the magnetic field 
evolution of the Andreev spectra, reproducing the main 
experimental features without any additional parameter 



(see Fig. EUc)). The theoretical analysis of the critical 
field reveals that for d > Ao, as in our case, both A and 
the spectral gap are finite up to the transition to the nor- 
mal state. This naturally explains why this transition is 
of first order and why the zero-bias conductance is not 
modified by the field. The existence of this first order 
transition in superconducting films was first discussed in 
the frame of the Ginzburg-Landau theory [2JJ| . 

In conclusion, we have presented a comprehensive 
experimental study of the transport through Al/Co 
nanocontacts. We have also introduced a model for the 
description of the Andreev reflection in S/F interfaces. 
While retaining the simplicity of BTK-type theories, our 
model includes the effect of a spin-dependent transmis- 
sion and allows the analysis of a great variety of realistic 
ingredients. We have shown that such a model consis- 
tently describes the whole set of measurements for ar- 
bitrary voltage, temperature and magnetic field, which 
demonstrates that the current polarization in ferromag- 
nets can be determined using Andreev physics. More- 
over, our data and analysis provide important input 
for first principle calculations of electron transmission 
through ferromagnetic interfaces. 
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